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Abstract 

<N 

In this work, we analyze a class of Yang-Mills models containing adjoint Higgs 
fields, with SU(N) symmetry spontaneously broken down to Z(N), showing they 
y—i contain center vortices, Y-junctions formed by them, and junctions where differ- 

ent center vortices are smoothly interpolated by monopole-like configurations. 
In the context of dual superconductors, these objects represent different states 
of the gluon field. Center vortices confine quarks to form normal hadron states. 
t-h The interpolating monopole, which in our model cannot exist as an isolated con- 

figuration, is identified with a confined valence gluon. A junction containing a 
monopole can bind quarks in a color nonsinglet state to form an overall neutral 
object, identified with a hybrid hadron. These states, formed by quarks bound to 
a valence gluon, are allowed by QCD, and current experimental collaborations 
are aimed at identifying them. Finally, considering the general version of the 
model, based on a compact simple gauge group G, the picture is completed with 
a heuristic discussion about why it would be natural using as G the dual of the 
chromoelectric gauge group G c , and external pointlike monopoles to represent 
the mesonic and baryonic Wilson loops. 
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1 Introduction 



Nowadays, the implementation of sophisticated lattice numerical techniques together 
with ever-expanding computing capacities is permitting to unveil the complex dynamics 
of the strong interactions. 

Using the tt, K and H masses as input parameters, the obtained light hadron spec- 
trum is in very good agreement with the experimental one PQ. The properties of the 
potential between heavy qq probes in pure SU(N) Yang-Mills (YM) theory were ob- 
tained for different representations of the external quarks. Among them, linearity and 
iV-ality at asymptotic distances, Casimir scaling at intermediate ones, and string-like 
behavior, have been observed [21 E] • Recent lattice calculations permitted to elucidate 
the type of interaction between three heavy quarks in pure SU (3) YM theory. The 
controversy about whether the 3q potential is A or Y-shaped at asymptotic distances 
was resolved in favor of the latter [4|-[8J. 

Regarding the properties of the gluon field, there is another central question in 
hadron spectroscopy, which is still open, concerning the existence of a color nonsinglet 
qq' pair neutralized in color by an excited gluon field, a state known as a hybrid meson. 
These qgq' states, where the quarks are bound to a valence gluon, as well as similar 
hybrid baryonic states, are allowed by QCD and have been searched for many years 
now. Among the new experiments, the glueX collaboration at the Jefferson Lab aims 
at mapping the spectrum of gluonic excitations using polarized photons. For a review 
of theoretical predictions for masses and decay modes, recent experimental evidence, 
and some lattice results, see [H]-[II] and references therein. 

On the theoretical side, motivated by the pioneering work of refs. [12j-[15j, Abelian 
Higgs models to describe the behavior of the qq potential [16] and the interaction 
among three quarks in normal hadrons, have been proposed in [171 [18] and |19j-|21j. 
respectively. Many efforts have also been devoted to understand dualities in supersym- 
metric theories containing non Abelian monopoles [22]-[27] and non Abelian vortices 
[28] 129] . for a review, see refs. [30l I3T] . In refs. [32]- [31], a non Abelian complex 
formed by an open vortex ending at a monopole/antimonopole pair, whose positions 
are kept fixed, has been obtained in supersymmetric models with hierarchical gauge 
symmetry breaking. In ref. [33], a Z(N) (2 + l)d model was constructed where not 
only quark/antiquark pairs are confined by domain walls, which in planar systems are 
stringlike, but also contains baglike baryonic states. Planar models that admit stable 
static domain wall junctions have been proposed in ref. [36J. 

In pure SU(N) YM theory, the iV-ality dependence of the qq potential points to 
center symmetry as a natural ingredient in the possible dual effective models. In 
particular, heavy quarks in the adjoint representation should not be confined [21 E], a 
property that would be naturally implemented by the triviality of vortices with center 
charge N. In ref. [37], the general discussion about the magnetic weights of non 
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Abelian monopoles [22] was adapted to classify center vortices in theories where the 
gauge group G is spontaneously broken down to its center (for a review, see [30J). An 
SU (N) model contaning stable center vortices, with at least iV adjoint Higgs fields, 
has been previously proposed in refs. [38j-|41j. 

In spite of all the above mentioned theoretical activity, a (3+l)d non Abelian model 
containing the different gluon field states needed to form normal and hybrid hadrons, 
is still lacking. In this article, we shall propose a first step towards this objective, 
concentrating on the underlying topological features. Although we will be mainly 
interested in SU(N) Yang-Mills-Higgs (YMH) dual models, most of our discussion will 
also apply to a general compact simple G with nontrivial center. 

In our model, the normal gluon field will correspond to center vortices, or Y- 
junctions exclusively formed by them, while the unconventional field, characteristic of 
either mesonic or baryonic hybrids, to different center vortices smoothly interpolated by 
monopole-like configurations, which represent the confined valence gluons. In addition, 
the quark probes will be represented by external pointlike monopole sources. 

Our work is motivated by a different but related aspect of gauge theories, namely, 
the search for the relevant configurations to compute averages in pure YM. In the lat- 
tice, some scenarios involving ensembles of magnetic defects exist, having their coun- 
terpart in the continuum. Monopoles have been introduced by means of a gauge fixing 
where a field transforming in the adjoint is diagonalized [12]. They have also been 
discussed in a gauge invariant way, in terms of a field decomposition that permits a 
simple description of monopoles as defects of a local direction n\ in color space [13] - 
[15] . Thin center vortices have been parametrized in [13 HB]- In ref. [1H], in the case 
of SU(2), we have written these singular objects as defects in a pair of local directions 
ri2,n3, which together with ni form a local color frame for the Lie algebra. We also 
showed this type of description extends to correlated monopoles and center vortices, 
configurations that in the lattice are known as chains. 

Now, the general discussion about the group weights and the topology for Z(N) 
vortices in dual superconductor YMH models [37] is similar to that for magnetic Z(N) 
thin center vortices in pure YM [UJ-[1H]. Therefore, a part of the analysis for YMH 
hybrid states should be similar to that for YM correlated vortices and monopoles. As 
we will see, in the YMH model defined by a gauge group G, the defects of a local 
frame for the Lie algebra q will provide the guiding centers for the smooth hybrid 
configurations. 

For this aim, the Higgs potential will be chosen such that the vacua are labeled by 
Ad{G), the adjoint representation of G. As a concrete example, we shall consider as 
many adjoint Higgs fields as the dimension of q, propose a flavor symmetry with global 
group Ad(G), and construct a potential whose terms only involve the product in the 
algebra. However, the general discussion could also be applied to other potentials with 
the same symmetry breaking pattern, and diagonal and off-diagonal sectors, setting 
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a scale for the monopoles and center vortices. This type of model has also points 
of contact with the SU(N) YMH model in [10], describing center vortices by means 
of adjoint Higgs fields along directions orthogonal with respect to a diagonal sector, 
although in that case the diagonal sector has been frozen in a way it does not contribute 
to the energy. The flavor symmetry also plays an important role in models containing 
non Abelian vortices [28], constructed with as many fundamental Higgs fields as the 
dimension of the fundamental representation of G. 

Regarding the two different aspects related with center vortices and monopoles, 
namely, their study in either the context of dual superconductor YMH models with 
gauge group G, or of pure YM theory with chromoelectric gauge group G e , the following 
remarks are in order. First, it is important to underline that, in the former case, a single 
configuration has a direct physical meaning, representing the chromoelectric confining 
string or junction acting between heavy quark probes. On the other hand, in the 
latter, a single configuration has no direct physical meaning, but represents possible 
virtual processes in 4D Euclidean spacetime, to be summed in physical objects such 
as the partition function. Second, while in YMH models the topological objects are 
smooth and essentially classical, in pure YM theory they are initially singular, but 
lattice results lead to consider the possibility they could become thick and relevant 
due to quantum effects [21 [3]. The hope is that in pure YM the effective action of 
magnetic topological configurations become small or negative [50]. In that case, they 
would proliferate and lead to a dual superconductor picture where similar topological 
objects, in a Higgs phase, represent the different states of the confining gluon field. 
For further considerations about this important point, see ^8] and refs. [5Tj-|53j. 

The article is organized as follows. In §[2j we present our model and analyze the 
possible vacua. In ^3] and §[4], we discuss the topology and charges of center vortices and 
junctions. Some basic properties of the Cartan decomposition of SU(N) are reviewed 
in £j5] In §[6j we show how to close the equations of motion for SU(2) and SU(3) center 
vortices in terms of a reduced set of profile functions, while in £j7] we discuss general 
properties of junctions and the associated non Abelian phases. In £|8j the natural choice 
of gauge group and quark sources in the YMH model is suggested. This is done in the 
framework of lattice inspired effective models, which also deal with topological defects, 
but on the side of the pure YM theory we are intending to describe. Finally, in £[9] we 
present our conclusions. 

2 The Yang- Mills- Higgs model 

As is well known, the SU(2) YMH model with a Higgs field in the adjoint, in the 
phase where the symmetry is spontaneously broken down to U(l), leads to the t' 
Hooft-Polyakov monopole. When a pair of adjoint Higgs fields is considered, the SSB 
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phase only preserves Z(2), and the model contains Z(2) vortices [38]- jH]. Motivated 
by the description of correlated singular Z(2) vortices and monopoles in SU(2) YM 
in terms of a local color frame m, 712, |19], it is expected that an SU(2) YMH 
model containing three adjoint Higgs fields will lead not only to smooth Z(2) vortices 
but also to smooth junctions formed by a pair of vortices attached to a monopole-like 
configuration. 

In this chapter we shall start by defining this type of model, generalized to the case 
of a connected compact simple Lie algebra q. The action is given by, 

S = Sym + Sfiiggs , (1) 
which contains the Yang-Mills term, 

Sym = -J d A x^ u -F^ , (2) 

F^ = F^-f={i/g)[D^D u ] , D„ = d ll -igA li , (3) 
where = ■ T, and the TVs, A = 1, . . . , d are hermitian generators of q, 

[T A ,T B ] = if ABC T c . (4) 

In the algebra, a positive definite metric exists, 

(X,Y)=Tr(Ad(X)Ad(Y)jT\, (5) 

where Ad(X) is a linear map of X e q into the adjoint representation, generated by 
d x d hermitian matrices M4, with elements M a \bc — —ifABC, satisfying, 

[M A , M B ] = if ABC M c , , Tr{M A M B )=5 AB . (6) 

With this normalization, we have, 

(T A ,T B ) = 5 A b , f abc Jdbc = $ AD ■ (7) 

The Higgs sector is given by a set of fields ipi = ipi ■ T, transforming in the adjoint, 

ijj^u^p^- 1 , ueG, (8) 

and the Yang-Mills-Higgs action is gauge invariant under the regular transformation, 

A^ -> Al = UA^U- 1 + -Ud^U- 1 , (9) 



1 If we were using a convention with antihermitian generators, in order to work with a positive 
definite metric, a minus sign should be included. 
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when accompanied by local transformations (|8|). That is, 

(10) 



D^i = d^i - ig[A^ V/] • (11) 

To construct a Higgs potential, defined on the fields ipj = if)j ■ T, the natural 
invariant terms involve the metric applied to pairs of Lie algebra elements. Note also 
that as we are using hermitian generators, the closed product in the algebra is in fact 
given by, 

& A ipj = -i[ipi,tpj] ■ (12) 
Then, up to quartic order, the terms we shall consider are, 

(ij>i,il>j) , (4>i, i>J A iPk) , (ipi A ipj, ipx A ipL) , (ipi,'ipj)(iliK,'tpL) ■ (13) 

At this point, an additional symmetry will be introduced, acting on the flavor, with 
global group Ad(G). For this aim, we note that the action of the gauge group in eq. 
8| can be represented using the components of ipi in the Lie basis, 



i>i = HaT a -» UtPjU- 1 = ^IbUTbIJ- 1 . (14) 
As the transformation of the TVs is a linear combination of generators, we can write, 

UTbU' 1 =T a Rab , (15) 

Ha -»- Rab $i\b , R = R{U) e Ad(G) , (16) 

where R is the d x d matrix that corresponds to the adjoint representation of the 
group. Then, to implement the additional symmetry, the index / will run from 1 
to d and, renaming / — > A, we shall require the action to be invariant under global 
transformations acting on the flavor index, 

ipA ^Hab^b , n e Ad(G) . (17) 
Up to quartic order, this restricts the possible terms in the potential to be Q 

u 2 

Vm ggs = c + — (ip A , ip A } + 

K A 

+ 3 f ABC ^Pa A i>c) + ^(^A tp B , fpA A ip B ) ■ (18) 



2 A (ipAiipB) 2 term could also be added, although the subsequent analysis would not be essentially 
modified. 
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Taking the first order variation of the potential with respect to ipA — > + ^Ai the 
critical points are obtained from, 



fi 2 ip A + k f ABC ipB A ipc + A ipB A (ip A A ips) = . (19) 
The possible solutions have the form, 

i) A = vST A S~ 1 , SeG, (20) 



with the condition 



fi 2 v + nv 2 + \v 3 = , (21) 



K / / K \ 2 jl 2 



that is, v — 0, or 

" = "«=-2A ± Vl2AJ -X- (22) 
Which configurations correspond to local or global minima depends on the choice of 
parameters. The associated Higgs potential is, 

Vffi gg s = c+^v 2 [(v- v ) 2 + b 2 } , (23) 

so in order to have an energy bounded from below, we must always have A > 0. We 
also note that changing in the potential k — > —n is equivalent to v — > —v, so we can 
start analyzing the k < case, and then extend the result to positive values. 

For control parameters in the range fi 2 > there is only one symmetric minimum 
at ipA = 0. When |^ > fi 2 > the global minimum is still at ipA = 0, and minima 
with higher energy appear at ipA = v c STaS -1 , with v c > computed using the positive 
root in eq. \22\. At /j 2 = l\, these minima also become ground states, degenerate 



with the symmetric one, and for /i < ~ their energy is lowered. In the SSB phase, 
the space of global minima is labeled by the adjoint representation of the group, 

M = {ip A = v c n A , n A = ST A S-\ SeG}, (25) 

where the n^'s form a local Lie basis, 

[n A ,n B ] = if A Bcn c , (n A ,n B ) = S A b ■ (26) 
In components with respect to the basis Ta, we can also write, 

M = {$ A = v c fi Al n A = R(S)e A , R e Ad(G)} , (27) 
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where eU is a d x 1 matrix, with its A-th element equal to 1, and the rest equal to zero. 

When the system chooses a given vacuum, say S — I, the continuous [/-symmetry 
of the potential is spontaneously broken to a discrete center symmetry ipA —> UiPaU^ 1 , 
U E Z(G), so all the components of the gauge field become massive. In addition, the 
continuous global symmetry 

*p A -> RabU^bU- 1 , R = R(U), (28) 

that involves global gauge and flavor transformations, remains unbroken. 

Note that unlike the regular gauge transformations R(U), when dealing with topo- 
logical objects, the extension of the x-dependent vacua to the whole space will be 
associated with a non Abelian phase R(S) containing defects. These will correspond 
to the guiding centers of vortices and monopoles. 

For k > 0, we have the same picture, with the minima at v c < corresponding to 



the negative root in eq. (22), and at k — 0, jj, < 0, there is a line where two types of 



degenerate ground states of the form {ip^ = v± ST^S 1 , v± = ±a/— /i 2 /A} coexist. 
2.1 Cartan decomposition of q 

In the next sections we shall discuss center vortices and junctions. For this aim, it will 
be useful considering the Cartan decomposition of the Lie algebra. 

A compact connected simple g can be decomposed in terms of diagonal generators 
T q , q — 1, . . . , r, which generate a Cartan subgroup H, and off-diagonal generators E a , 
or root vectors of the adjoint representation, labeled by the roots a = (at, . . . ,a r ), 
satisfying, 

[T q ,T p ] = , [T q ,E a ]=a q E a , [E a , E_ a ] = a q T q . (29) 

[E a , E y ] = N ai E a+1 , a + 7 ^0, (30) 

where N ai = 0, if a + 7 is not a root. The off-diagonal generators can be grouped 

in pairs E a ,E_ Q , where a are the positive roots. The hermitian generators, can be 
identified with, 

{TJ , {TJ = {T a , T a } , a > , (31) 

T a = -j={E a + E- a ) , = -j=-(E a — E^ a ) , (32) 
where the following triplets generate so (3) subalgebras, 

« • T , -f= T a , — L= T a . (33) 



a Va 2 Va 2 

Note that for the set of all hermitian generators we are using upper-case letters, 
Ta, reserving lower-case letters, in the beginning of the alphabet, for the off-diagonal 
sector, T a . 
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The weights of a given matrix realization of the Lie algebra will also play an impor- 
tant role. A weight w is defined by the eigenvalues of diagonal generators corresponding 



to one common eigenvector. In particular, the second equation in (29) says that the 
roots are the weights of the adjoint representation, which acts via commutators. We 
shall also refer to the eigenvectors as weight vectors. 

2.2 Local minima 

At the classical level, we shall be interested in looking for topologically stable (dual) 
magnetic static solutions to the equations of motion, whose total energy is given by, 

E = Jd 3 x (jiF^Fij) + i(A^A, AV»a> + V mggs (ip A )j . (34) 

We are particularly interested in the SSB phase, where the global minima are given 
by Ad(G). In this phase, let us intially analyze the local ground states, that is, the 
general conditions necessary to have, in a given asymptotic region, energy densities 
tending to zero. When discussing static configurations, we shall partially fix the gauge 
by the condition A = 0. More generally, we could also be interested in analyzing the 
path integral quantization of models containing a Yang-Mills term in 4D Euclidean 
spacetime, where relevant configurations must have finite action. In that case, the 
asymptotic conditions on the gauge field A^ are simply obtained by replacing i — > // 
in the equations below. 



In the asymptotic region, not only the condition (25) but also DiipA = must be 
imposed, 

DiUA = , ua = STaS^ 1 . (35) 

The most general A\ can be locally written as A{ |i oc = Af , a gauge transformation 
satisfying Di(A) Ta = 0, that is, 

[Ai,T A ]=0. (36) 
As this is valid for every A = 1, . . . , d, and the algebra is simple, this implies, 

A = , Aiho^-SdiS- 1 , (37) 

9 

leading to F^- = on that region. 

We would like to underline that in general, even in the asymptotic region, eq. (37) 
is only valid locally, as S could change by a center element when we go around a 
loop. How to extend this expression to the whole asymptotic region depends on the 
global structure of G used to write the model. However, as all the fields are in the 
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adjoint, any structure represents the same model. One possibility, when using a general 
representation S, is to write the asymptotic configuration in the global form, 



Ai = -SdiS' 1 - Ii , (38) 
9 

where Ii is concentrated on the surfaces where S changes by a center element, and 
is designed to compensate the derivative of the discontinuities in the first term. This 
type of approach has been introduced in [47J [48] in the context of pure YM theories. 
This is simplified when using the adjoint representation of the group, which does not 
"see" the center of the group, 

AfM A = -RdiR- 1 , R = R(S) , (39) 
9 

as R will always be single-valued when we go around a loop. Another convenient 
manner to deal with this problem, in any representation, is by writting the asymptotic 
gauge field in terms of the single- valued quantities n A = STaS' 1 [4"9"] . [SI]- [S3]- In this 



respect, using that in the asymptotic region Ai is defined by eq. (35 ), with A = 1, . . . , d, 
we have, 

[n A , D inA ] = , (40) 

and using the property, 

n A A (X A n A ) = (X, n B ) n A A (n B An A ) = I, (41) 

we get, 

A i = -(l/g)n A Ad l n A . (42) 
Alternatively, expanding d^A = (nB,din A ) n^, we can also write, 

A = -Cfn A , Cf = -(l/g)f A Bc(n B ,d i n c ) (43) 



(see also ref. 

This can be extended to the whole space, so that the gauge field Ai, whose general 
form is locally given by, 

Al^^SAiS-'+'-SdiS- 1 , (44) 
9 

can be globally written in the forms, 

AfM A = R (AfM A )R~ 1 + l -RdiR- 1 , R = R(S) , (45) 
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A t = (A?-Cf)n A , Ai = Af T A , (46) 

with the condition, 

Af , (47) 

when approaching the asymptotic region. At the guiding centers of the topological 
objects, some frame components n A will contain defects and, as a consequence, some 
of the frame dependent fields Cf will be singular. For this reason, we shall need to 
impose additional boundary conditions on the gauge and Higgs fields. In fact, we will 
see in £3.1 ^6] and ^7] that among the quantities Cf, some diagonal components play a 



similar role to d^ip, when analyzing a Nielsen-Olesen vortex, while some combinations 
of the off-diagonal components, a similar role to f x <9jf , when analyzing a 't Hooft- 
Polyakov monopole 



To compute the field strength, we can use eq. (45) and recall that the commutator 



[di, dj] is nonzero when applied on singular quantities, to obtain, 

F$M A = R(F£(A)M A )R- 1 + (i/g)R[d i ,d j ]R- 1 , (48) 
where the second term can also be rewritten as —(i/g)R(R~ 1 [di, dj]R)R~ 1 . In addition, 



from eq. (43), we have 



CfM A = -R~ 1 d l R, , {i/g)R~ 1 [d i ,d j ]R = F l A {C)M A . (49) 
9 



Then, putting these pieces together, we get, 

F ij = F^T A = Gf j n A , G? 3 = F$(A) - F$(C) , (50) 

where F^{C) is concentrated at the frame defects [5^] . 

F*(C) = (i/g)tv(M A R- 1 [d i ,d j ]R) 

= d t Cf-d 3 C? + gf ABC C*Cf. (51) 

2.3 The Yang-Mills-Higgs equations 



The minimization of the energy (34) gives, 

DjFij = ig tyu, D*I> A ] , (52) 

DiDiip A = n 2 ?p A + k f ABC ipB A ipc + A ipB A (ip A A ip B ) . (53) 

For a given map S, defining a local frame n A containing defects, we can look for 
solutions with the gauge fields given by eq. ((461 and the Higgs fields by, 



ip A = S{h AB T B )S-\ (54) 
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where h AB are real symmetric smooth profile functions. 



Using eq. (44), or its global form (45), at every point but on the zero measure 



region of frame defects, we get, 

VfF?{A)=gf AB ch DB D° E h DE , Vf B = 5 AB d i -gf ABC A? , (55) 

T>f C T>f D hAD = fJ 2 h AB + K f BCD f A.EF ^EC^FD 

+ ^ jBCofoEF hcch A EhGF ■ (56) 

In the asymptotic region, besides Ai —> 0, we have to impose, 

h AB -> Vc $AB , (57) 
and, at the guiding centers, h AB and Ai must satisfy boundary conditions such that Ai 



and ip A , in eqs. (46) and (54), be smooth. These requirements will be further modified 
by the addition of external quarks sources. 

In £[6j we shall give a center vortex ansatz such that the equations of motion can 
be written in terms of a single profile function for the gauge field components, and a 
reduced set of Higgs profiles. Before doing this, general model independent properties 
about the topology, weights and charges of vortices and junctions are developed in the 
following sections. 



3 The topology of center vortices 

In the absence of external quark sources, the relevant global structure of the model is 
that of the adjoint. It can be given by Ad{G) = G/Z(G), where Z(G), the center of G, 
is the kernel of the homomorphism G — > Ad{G). This is a discrete subgroup formed 
by those elements S G G verifying ST A S~ 1 — T A , for A = 1, . . . , d. 

In particular, to compute the fundamental homotopy group of the adjoint, it is 
useful considering the universal covering group G, uniquely defined as the simply con- 
nected group whose Lie algebra is jj, and apply the exact homotopy sequence, 

= IIi(G) -> U^G/ZiG)) -)• U (Z(G)) -> U (G) = , (58) 

to conclude that the fundamental homotopy group of the space of vacua Ai in eqs. 



(|25]), (|27j) is, 

Ux(M) = U 1 (Ad(G))=Z(G) . (59) 

For example, when G = SU(N), Z(G) = Z(N). When we go along an open path 
S(t) E SU(N), t G [0, 1], that changes by a center element e l27r ^ N , a nontrivial closed 
path in the adjoint representation of SU(N) is obtained. If this path is composed N 
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times we get a trivial map in the adjoint, as it is associated with a closed path in the 
fundamental representation of SU(N), whose first homotopy group is trivial. This is 



the meaning of Tl^AdiSU (N))) = Z(N) in eq. ((59]). 

For a straight center vortex, the asymptotic region is given by R 3 — T 2 , the comple- 
ment of a cylinder T 2 containing the vortex. The defining property of a center vortex 
is that for any closed path contained in the asymptotic region, linking the vortex core 
T 2 , the Wilson loop, 

W[A] = (l/d)trPexp (ig j> dx^AA , (60) 

computed using the group elements in G, gives a center element in Z(G) (P stands 
for path ordering). Note that this is a gauge invariant definition, as for regular gauge 
transformations W / [v4 c/ ] = We would also like to underline that it is only in 

the context of pure Yang-Mills theories, where the gauge field represents the gluons 
and the path a quark/antiquark pair, that the Wilson loop is a direct measure of 
the interquark potential. In the dual YMH model, in order to represent the quarks, 
appropriate external pointlike monopole sources will have to be introduced, see the 
discussion in £j8j 

Now, as in the asymptotic region the Wilson loop assumes a discrete value (the 
group is simple) it cannot be changed by continuous deformations of the smooth local- 
ized ansatz, as long as the loop be always chosen contained in the asymptotic region of 
local minima. Then, when Ui(Ad(G)) is nontrivial, there will be static center vortex 
configurations with finite energy per unit length of the vortex. 



3.1 Center vortex weights 

Let us consider a vortex characterized by a center element 3 G Z(G). This means that 
the asymptotic gauge field must be such that, 

W[A] = 3 , (61) 

where C is a loop completely contained in the asymptotic region, linking the vortex 
core. This is the case when using in eqs. (37), (39), 

S = ? R = ^0-M^ (g 2 ) 

with the possible /3's satisfying, 



(63) 
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thus leading to a single- valued R transformation in Ad(G), 



R(2n) = e"* p - M = I dxd . (64) 

Here, (p is the polar angle with respect to the center vortex axis, and we adhered to the 
convention that in dot products involving a weight, (3 is considered as . . . , /3 r , 0, . . . , 0). 
The possible center vortex weights have been discussed in refs. jlOl |37] , following 



the ideas introduced in [22] to characterize non Abelian monopoles. Using eq. (64) on 
the weight vectors of the adjoint representation Ad(G), we get, 



e WQ = 1, (65) 

where a is a weight of the adjoint, that is, a root of $j. Then, for every a, the possible 
weights must satisfy the condition, 



whose solution is given by 



(3-aeZ . (66) 



i G A(G V ) • (67) 



That is, when the Higgs fields transform in Ad(G) = G/Z(G), the possible /3's are 
associated with the weights of the universal covering of the dual group G v , whose Lie 
algebra is defined by the dual root system, 

a< = * 5 • (68) 

For simple Lie algebras, N is a number, or rescaling parameter, that only depends on 
the general structure of q. In the lattice A(G V ), there are many types of center vortices 
labeled by 3, although this value does not fix uniquely the asymptotic behavior. A pair 
of asymptotic non Abelian phases defined by /3 and /3' are topologically equivalent if 
the map 

RR >-l = e i<P0-0')-M^ (6Q) 

can be continuously deformed to the trivial one. This is the situation when the asso- 



ciated path in the simply connected G is closed. Then, substituting M in eq. (69) by 
generators of the universal covering G, and acting on the possible weight vectors, two 
vortices are topologically equivalent if, 



e 
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where w is any weight of the universal covering G. This means that the closed path 
in Ui(Ad(G)) will also be closed in G that, being simply connected, will permit the 



continuous deformation of the composition (69) to a point. The solution to eq. (70), 
i.e., ((3 — (3') ■ w G Z is given by, 

i 0- $) e A(Ad(G v )) , (71) 

the lattice of the dual root system [22]. 

In other words, a center vortex defined by 3 can be associated with equivalence 
classes of asymptotic behaviors, 

Z(G) = A(&)/A(Ad{G v )) . (72) 

The weights of a given irreducible representation only belong to one equivalence class 



An important point we would like to emphazise is that, for a given 3, we are 
interested in determining the smooth configuration that minimizes the energy. If we 
start with an asymptotic behavior labeled by /3', the search for this configuration should 
not only contemplate the direct extension of the phase inside the core, keeping it along 
a Cartan direction, but also other possibilities. The behavior S = e tifl3 '' T on a circle 
contained in the asymptoti c reg ion could also be continuously changed to S = e llpl3 ' T , 



with ^ /3' satisfying eq. (71), when approaching the vortex axis. In that case, the 
extended phase S(x) would use the whole non Abelian group. 

For example, an asymptotic behavior labeled by a nonzero weight in A(Ad(G v )) 
can be continuously converted into the behavior S ~ / at the z-axis. Then, in this 
case, there are in fact no defects, and the minimization process will simply return 
a trivial result, corresponding to a pure gauge transformation of the configuration 
Ap = 0, ipA = v c Ta- On the other hand, when the asymptotic weight is not in 
A(Ad(G v )), any phase extension inside the vortex will necessarily lead to a local frame 
Ha = STaS^ 1 containing defects. For those /3 values characterizing the behavior of the 
stable solutions at the origin, it will be possible to perform a gauge transformation such 
that the phase is everywhere along the Cartan sector. In that case, it will be useful 
to compute the Abelian charge of the center vortex. For the minimum charge center 
vortices in SU(N), these values are expected to be the weights of the fundamental 
representation [37J. 

3.2 Charges of center vortices 

When the smooth configuration minimizing the energy contains the Abelian phase 



factor (62) at every point x 6 Br, the local frame ua, the frame dependent field Cf in 
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Figure 1: A Z(2) vortex. The red arrows represent the diagonal direction, while the 
green and blue arrows, the off-diagonal local ones. 



eq. (43), and F-(C) can be easily computed. For example, using eqs. (49), (51), the 



nonzero components are the diagonal ones, A = g, 



C\ = --d^^ , F? j (C) = --[d i ,d j ] l pf3\«. (73) 

In that case, in order for the gauge field and field strength be well defined smooth 
functions, the condition, 

A\ -> CI , (74) 

must be satisfied at the center vortex axis. In addition, the off-diagonal components 
that rotate when we go around the z-axis, 

ST a S~ l , n a = STaS -1 , (3 ■ a ^ , 



(75) 



contain a defect at the origin, so those profile functions multiplying them should tend 
to zero there. In Fig. 1, we show the local frame for a unit charge center vortex, in 
the case of SU{2). 

The center vortex charges have been obtained in [37] . Considering the dual tensor 
components along the local direction ha, 



rA — c nA 



jk ) 



(76) 



and organizing the fluxes Ql = J dSi Q\ as an r-tuple, Q v = (Ql, . . . , Ql), we get, 



Qv = —0 

9 



(77) 



These charges are invariant under regular gauge transformations (see the discussion in 
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4 The topology of junctions 



As the second homotopy group of a compact group is trivial, we have 

U 2 (M) = U 2 (Ad(G)) = , 

that is, there are no isolated monopole configurations. However, we will see that center 
vortices can concatenate some monopole-like objects to form junctions. These will be 
nontrivial topological configurations along a local Cartan sector of the Lie algebra, 
necessary to interpolate the center vortices originated by SSB. In this sense, it is more 
appropriate referring to them as monopole-like junctions, as their existence depends 
on the attached vortices. 

While a single center vortex is associated with smooth extended local Cartan di- 
rections h q {x) = R(x) e q , this will not be the situation in a monopole-like junction. In 
order to understand them, let us first discuss the su(2) case, and then generalize to q. 
In the framework of pure YM theories, we have seen in refs. jl9] that when using field 
decompositions involving a local color frame ha, A = 1,2,3, chains formed by corre- 
lated monopoles and thin center vortices can be represented by thin junctions formed 
by correlated defects. In su(2) there is one Cartan generator T 1; and given a matrix 
R E Ad(SU(2)) = SU(2)/Z(2) = SO(3), an associated local direction h x = Re x . On 
the other hand, given a unit vector hi, R is defined up to right multiplication by a 
matrix leaving e\ invariant. 

Now, let us suppose we give a well defined configuration hi(x), with x on a surface 
S 2 around a given point. There are two possibilities, the first one is when this config- 
uration can be generated by means of a well defined R(x), h\{x) = R(x)ei, x E S 2 . 
In that n 2 (S'0(3)) = 0, R(x) can be continuously deformed to the trivial map 

R = I, and in this process the configuration h\{x) will be continuously deformed to the 
trivial one h\(x) = e%. However, there is the possibility that a well defined configura- 
tion hi(x) on S 2 be generated by an ill-defined R(x). In such a case, we could not use 
the argument that n 2 (SO(3)) is trivial, to imply a trivial hi(x), as R(x) would not be 
defined on S 2 . Take for instance the configuration shown in Fig. 2, where, we have a 
hedgehog form for h\(x) (red arrows), well defined on S 2 . This is possible as there are 
two points on the sphere where the off-diagonal directions (green and blue arrows) are 
ill-defined. When we go around a small loop around the north (south) pole, h a = Re a , 
a = 2, 3 rotate once, so they are ill-defined at the poles. These rotations have the same 
orientation with respect to the outward normal, but different orientations with respect 
to z. 

The presence of singularities can be understood in terms of the hairy ball theorem. 
When h\ is along the radial direction, hi and n.3 are tangent to S 2 , so these fields 
must be ill-defined at a pair of points, namely, the guiding centers of the vortices. The 
non Abelian phase we have described then corresponds to a pair of center vortices, 
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Figure 2: A point-like monopole at the border of a pair of open Z(2) vortices. 



attached to a nontrivial configuration in the coset SU(2)/U(1) = S 2 = {hi}. Of 
course, to extend this type of local frame to a well defined field configuration on the 
three dimensional space, we still have to impose appropriate boundary conditions. On 
the z-axis, the profile functions multiplying fi2 and should tend to zero. Then, 
when looking at the regular fields if) a on the sphere S 2 , they will only form a basis 
at points other than the guiding centers of the vortices, this is the reason why these 
configurations are not governed by n 2 (5'0(3)) =0. In addition, the profile functions 
accompanying h\ should tend to zero at the origin, where the hedgehog h\ is ill-defined. 



For general G, to construct the junction, we can start by giving a well defined 
configuration for the local Cartan directions n 9 , on a two-dimensional surface S 2 around 
a point. Again, when representing this configuration by n q = STgS -1 , or in components 



with respect to Ta (see eq. (27)) 



h g (x) = R{x) e g , x e S 2 , Re Ad(G) , (78) 

R(x) could be a well defined map on S 2 , and in that ) would be trivial, or 

it could be ill-defined. That is, on some disks D c G S 2 of small radius e, centered at 
points c = ci, C2, • • ., S(x) and R(x) could behave as, 

S c (x) ~ S c e iXc & f , R c {x) ~ TZ C e lXc ^ (79) 

(no summation over repeated c-type indices). Here, the first factors are smooth, while 
the angles Xc parametrize the borders dD c , oriented in the positive sense with respect 
to the outward normal to S 2 . The /3 c 's are in A(G V ), the set of possible weights leading 



to a single- valued R c (x) when we go around these borders (see { 3.1 ). With the behavior 
(79), the smoothness of n q (x) is guaranteed as the singular phases have no effect on 



Tlq at C , 



lim n q = S c T q S c 1 , c = Ci, c 2 , . . . . (80) 
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On the other hand, some off-diagonal components n a (x) will be ill-defined at the points 
c, which correspond to the guiding centers of the vortices attached to the monopole. 

An important global constraint on the possible field configurations is obtained by 
noting that, once the singular points are determined, the composition 



r = 71 o 0D Cl 07/0730 0D C2 o 7 2 1 



o 



is a closed path that can be deformed to a point in S 2 — (D Cl U D C2 U . . .), a region free 
of singularities. Here, 71, 72, • • • , are open paths on S 2 joining a common point Xq with 
points at dD cl ,dD C2 , . . ., respectively. Let us suppose the above mentioned region is 
contained in the asymptotic region, that is, the field is pure gauge there. Then, the 
Wilson loop computed in G, for any closed path contained in S 2 — (D Cl U D C2 U . . .), 
must give a center element in Z(G). Now, the center of a simple group is discrete, so 
that the value on T, 

W r [A] = (l/d)tr [ e *4' f e *4' f .. .] , ( 8 1) 

cannot change as T is deformed to a point, using paths contained on a region free of 
singularities. That is, we obtain the constraint, 



J2tt/3 c1 -T J2ir/3 C2 -T 



\2) 



^(i+i + ...)6A(MG v )). (83) 

This is a matching condition between the center vortices and the attached monopole. 
In other words, the composition of the loops e* Xc/3c ' M , c = ci, C2, ■ • • , is an element in 
Ad(H) that must be closed in G and, therefore, topologically trivial in Ad(G). 

The use of exact homotopy sequences is the appropriate framework to settle the 
discussion above, and obtain the conditions to have a topologically nontrivial vortex- 
monopole junction. These sequences have also been used to analyze non Abelian 
monopoles and vortices in systems with hierarchical gauge symmetry breaking [30] . 
The essential ingredient in our construction is that at any given point, the n q 's define 
R G Ad(G) up to right multiplication by an element in Ad(H). Then, the set of rz g 's 
can be identified with the quotient space Ad(G)) / Ad(H) . 

As is well known, for relative homotopy classes associated with a pair of manifolds 
X, Y there is an exact sequence: 

u 2 (Y) -> n 2 (Y, x) -)■ riipf ) -> iL(y) , (84) 

such that the image of a given map in the sequence is the kernel of the next map. This 
theorem can be applied to Y — Ad (G) and X = Ad (H), to obtain, 

= U 2 (Ad(G)) -> n 2 {^^yj -> III (Ad (H)) -> n x (Ad(G)) = Z(G) . (85) 
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The image of the first set in the sequence can only be a trivial map in the quotient, 
which coincides with the kernel of the second map. In other words, different equivalence 
classes in IL^Y/X) must correspond to different classes in Ux(X). If the last set in 
the sequence were trivial, then the kernel of the third map, and therefore the image of 
the second, would be the whole IIi(X). In that case U.2(Y/X) would be homomorphic 
to I1 1 (X). However this is not the case, there are classes in I1 1 (X) that cannot be 
attained by the second map. The only classes attained are those mapped into the 
trivial class in II 1 (y) = Z(G). 

The S 1 manifold in Ili(y) precisely corresponds to dD Cl o dD C2 o . . ., the closed 
path obtained by composing the closed paths around the center vortex defects. In 
this regard, note that in the simply connected case, the construction of the sequence 
starts by showing that any map G(s,t), from a square into G, is equivalent to a map 
where three of the sides are mapped to the identity J, while the fourth is mapped 
into H. Identifying the sides adjacent to the fourth, we get a map from D — D±, a 
disk with a hole at D\ C D, with the external border dD mapped into J, and the 
internal border 8D\ into H. As the points in 8D\ and dD can be identified in G/H, 
an element in U 2 (G/H) becomes correlated with an element in Ui(H). Note that in 
such configuration, the loop in Il 1 (if) must be trivial. When we continuously move 
from 8D\ to dD, it must be continuously deformed to the identity map. As G is simply 
connected, this condition is automatically met. 

However, when Y is not simply connected, for a loop contained in H at dD±, 
that is nontrivial in G, this construction is not possible. For Y = Ad(G), the best 
we can do is considering a disk D with a number of holes at Di,D 2 , . . ., where the 
external border dD is mapped into /, and on the borders dD\,dD2, ... we have maps 
into Ad(H), with the condition that their composition be trivial in Ad(G). This is 
the previously discussed junction composed by center vortices, with trivial total Z(G) 
charge, correlated with a monopole-like configuration in the coset Ad(G)/Ad(H) = 

When restricting to those configurations satisfying the matching condition, we have 
U 2 (Ad(G)/Ad(H)) = U.iAdmU . (86) 
Therefore, inequivalent elements in Ui(Ad(H)) correspond to inequivalent elements in 



Tl2(Ad(G) I Ad(H)) . Then, among the possible charges in eq. (83), the condition for 
having a nontrivial monopole-like junction is 

4 + (3 C2 + . . . + , (87) 

while configurations with /3 Cl + (3 C2 + . . . = are exclusively formed by center vortices. 
For example, in SU(2), a single center vortex (Fig. 1) and a pair of center vortices 
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attached to a monopole (Fig. 2), that correspond to zero and two turns of the off- 
diagonal components when we go around dDi odD 2 , cannot be continuously deformed 
one into the other. 



4.1 Charges of junctions 

For junctions containing a monopole, it is natural defining a monopole asymptotic 
region R 3 — S 3 , the complement of a ball, where the quotient does not contribute to 
the energy, that is, the fields are stabilized so they satisfy, 

ip q {x)=v c n q , Din q = , n q = ST q S~ 1 . (88) 

The most general Ai can be locally written as Af , which in the asymptotic monopole 
region must satisfy Di(A)T q = 0, that is, 



[Ai,T q ] = . 



(89) 



As this is valid for every q — 1, . . . , r, and the T q s form a maximal set of independent 
elements in the algebra, commuting among them, Ai must be a combination of the 



Ai = A\T q , A i= A\ STgS- 1 + -SdiS' 1 , 



(90) 



or defining Af = Af — Cf, and using eq. (46), this implies, 



Ai -> A\n q - C«n a 



(91) 



Then, using eq. (50), in R 3 — S 3 the local diagonal and off-diagonal components of the 
field strength are, 



= (diA'-d^-gf^CtC) 
G% = -F«(C). 



(92) 



From eq. (51 ), around a given vortex of the junction, located at c = ci, C2, • • ., we have 
(cf. eq. m), 



F A (C) = -tY{M A e- ix °^ M n- x [di,dMn c j x °^ M ) 

9 

= -i t r(M A /3 c -M)[^,^]xc, 



(93) 
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That is, in R 3 — S 3 the field strength is, 



G % n g , 



KMJ - djAf) + II 



-gr b CtC" . 



(94) 
(95) 



To put in evidence the constraints implied by the exact homotopy sequence (58), 



let us compute the charge of the monopole along the diagonal local direction n q . Using 



eqs. (|76|), (|95|), it is given by, 



dSj hi 



hi 



, TTl 



(96) 



where the flux is computed through a closed surface contained in the asymptotic 
monopole region. On the other hand, taking into account that a part of this re- 
gion is also an asymptotic region of local minima, where Fj k = 0, we can also compute 
the total flux as the sum of the fluxes through the disks D c , chosen so that the dD c 's 



encircle the vortex cores. That is, we can write, 



Q q m = J2j Mi g " = E J dS > M A l + h D 



(97) 



Now, taking into account eq. (93), the flux for each vortex is given by, 



l 



dSi - e ijk Ff k (C) 



- tr (M^ c ■ M) J dS, e ijk dfrx* = jPd 9 



so the possible monopole charges, organized as Q r 

2tt 



(Qln 



Qln), become 



Q, 



g 



(98) 



(99) 



They correspond to (2ir/g) 2N times a nonzero element of the dual root lattice A(Ad(G v )) 



(cf. eq. (83)). Note that among the possible monopoles, those with charge given by 
the roots of the dual algebra $j v are expected to be elementary, as these roots generate 
the lattice. 

Note also that a regular gauge transformation corresponds to R — » 7ZR, where 1Z 
is regular and R defines the frame defects. Then, as the singular term FMC) is not 



altered by this operation (cf. eq. (51)), the center vortex and monopole charge, along 



the local diagonal directions, are in fact invariant under these transformations. 
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5 SU(N) 



In this section, we review some basic properties of the Cartan decomposition, and the 
associated roots and weights of su(N). We shall follow ref. [55], but with a modified 
convention for the generators, so as to be consistent with the normalization we choosed 
in eq. 

5.1 Cartan decomposition of su(N) 

The rank of su(N) is r = N— 1, its dimension is d = N 2 — 1, and there are N(N — l)/2 
pairs of root vectors E a , E_ a . In the fundamental representation, the diagonal of T q 
can be given by, 

-^ W (1,...,1,- S ,0,...,0), (100) 

where the initial q elements are equal to 1. In this case, it is convenient to define a 
weight as positive if the last nonvanishing component is positive, so that the weights 
of the fundamental matrix representation, w q , q = 1, . . . , N, given by, 

ti q = (Ti\ qq ,T 2 \ qq , . . . ,T N _i\ qq ) , wi + . . . + w N = , (101) 

become ordered according to W\ > w 2 > ■ ■ ■ > Wn > and satisfy, 



N — 1 1 

*V^ = -^- > *V™p = -7^ . I^P- ( 102 ) 



The positive roots are, 

a qp = w q - w p , (103) 

with q < p; they satisfy a 2 = a ■ a = 1/N. This is a set of N(N — l)/2 elements 
that can be ordered in the form a\ > a.% > . . ., using the above mentioned notion of 
positivity. With our conventions, the sum over all the roots (positive and negative) is, 

^2a\ q a\ p = S qp . (104) 

a 

In the fundamental representation, the root vector E a is an iV x N matrix that, 
for a = a qp , only has a nontrivial element at position qp, whose value is l/y/2N. 

For example, in SU(2), the weights are one component. Those of the fundamental 
representation are uj\ = ^4=, o; 2 = — and there is a positive root a,\ = W\ — w 2 = 

l/\/2. The relation with the usual generators r A /2, where r A are the Pauli matrices, 
is as follows, 

Tl = 27I ' Tai = 271 ' T " 1 = 27f- (105) 
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In SU(3), the weights of the fundamental representation are, 

iBi = (1/2^3,1/6) , w 2 = (-1/2^3,1/6) , w 3 = (0,-1/3) . (106) 
Then, the positive roots are, 

3i3 = (l/2V3,l/2) , a 23 = (-1/2V3, 1/2) , 5 12 = (1/V3, 0) , (107) 
whose ordered form (c?i > c? 2 > ^3) is given by, 

«1 = «13 , "2 = "23 , "3 = "12 • (108) 

The relation with the generators A /2, where A A are the Gell-mann matrices, is given 
by, 

Ti = , T 2 = , (109) 

T ^ 4 7^ _ ^6 rp _ Al 



2V3 ' 2 2^3 ' 2^3' 

Ql 2v/3 ' a2 2^3 ' ° 3 2^3' 1 J 

Finally, a part of the 577(3) Lie algebra can be written as, 

[7q, Tq,] icx\^Ta , [Xjj,, Tg.] iot\qT a , [T^,,Tq,] icylqXjj , (111) 

corresponding to a structure constant = a\ q . The indices a and a refer to the 
same root a, but associated with the generators T a and T & , respectively. The remaining 
commutators, with a 7^ 7 are, 

[T Q , T 7 ] = if ai $T$ , [T a , T 7 ] = ifa^Ts , [Ta, ^7] = ifajgTs , (H2) 
with the nontrivial structure constants given by, 

/«i«2«3 faia2d3 /01O203 /cf3Cf2<5i 2\/5 ' (^^) 

and cyclic permutations. Here, the indices ai, a 2 and 03 refer to the roots given in eq. 
p7b. 
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5.2 Weights of Z(N) vortices 

For su(N), N is given by (d — r)/r = N, and a 2 = N, so that a v = a, and the possible 
/3's are (see [22], [30J and references therein), 

P = 2Nw , (114) 
where w are the weights of any matrix representation of SU(N). In particular, recall- 



ing that the weights of the fundamental matrix representation satisfy eq. (102), and 



the roots are given by their differences (cf. eq. (103)), the desired single-valuedness 



condition in eq. (65) can be verified for 



P q = 2Nw q , 



f3 q ■ a 



(3 q ■ a pp/ = (p,p ^ q) . 



(115) 
(116) 



As noted in ref. [37], fundamental and antifundamental weights correspond to 



the minimum charges. Applying SV,(27r) (cf. eq. (62)) on any weight vector of the 
fundamental representation, it is obtained, 



A2n B„-w v 



(117) 



27r(AT-l) 
N 



■W ( if ? 



— * n 

and as 2n (3 q ■ w p corresponds to either — (if q ^ p) or 
p), a topological charge z = — 1 is implied. The generators of the antifundamental 
representation —T* A have associated weights — w q , thus leading to a topological charge 
z = l. 



6 Center vortex ansatz 

For a straight center vortex with fundamental weight (3 q , let us consider the extended 

S = e iv0 «- f . (118) 



map (cf. eq. (62)) 



From eq. (29), we have, 

q -T,E a ] = (a-P q )E a , SEaS- 1 = e« 5 ^E a . (119) 
Therefore, for the extended local frame, besides the trivial diagonal components, 

n q {x) = ST.S- 1 = T q , (120) 

we get 

n a (x) = STaS^ 1 = cos(a • (3 q )<pT a — sin(a • (3 q )<pTa , 

ria(x) = ST a S~ l = sm(a ■ (3 q )(pT a + cos(a ■ (3 q )(pTa , (121) 
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and using eq. (116), for the positive roots of the form a = a pp >, p,p' ^ q, we have 
n a (x) = T a , n a (x) = T a . On the other hand, for those roots where the weight w q 
participates, the pairs n a (x), n a (x) rotate once when we go around the center vortex, 
so these frame components have a defect at the z-axis, where the guiding centers of 
the vortex are located. The cases a = a qp > (q < p') and a = a pq (p < q) have opposite 
senses of rotation. Then, we shall propose the following ansatz, 

ip q = h qp T p , ip a = h^ST^S' 1 , ipa = h^ST^S' 1 , (122) 
A = -a{p)d i i P P q ■ f . (123) 

which corresponds to, 

A = -(a-l)%>/VT . (124) 
From eq. ( 121[ ), it is natural to assume h a = h^, and use for h a only two profile 



functions h(p), h (p), depending on whether the components n a ,n a rotate or not. 
When p — > oo, in order to have a configuration of local minima, we require, 

h qp ->v c 5 qp h(p)->v c , h (p)->v c , a(p) ->■ 1 , (125) 

while the regularity conditions when p — > are, 

h{p) -> , a(p) ->■ . (126) 



Using this ansatz in eqs. (55), (56), and after putting h a = h a , the equations for 



the gauge field and diagonal profiles become 
1 fa' 



(?L- a ' ,y j (3 q -f = 2g(a- q ){l - a)h 2 a a-T, (127) 

(d 2 h qp )T p = p 2 h qp T p + 2(k a\ q + A a\ q >h qq/ )h 2 a a\ p T p , (128) 

with summation over repeated indices, q', p, and a. The sum over a is only done over 
the positive roots. From the field equation for if> a , we obtain, 

(d 2 h (a) - (1 - af{a- P q ) 2 (l/p 2 )h (a) )T ia) = 

= (p 2 h( a ) + 2Kh( a) h qql a\ q a\ q > + A hf a) a 2 )T( a ) 

+A h^h^Tj A (T( a ) A T 7 ) + A h 2 h^ a )T^ A (T (a ) A T^) 

+2nh 7 hs f(a)^sf^s(a)T{a) + A a\ q >a\ q "h( a )h qq >h qq >'T( a ) . 

(129) 

Here, the sum over 7 is done with 7 ^a, and there is no summation over a. As h a = h a , 
after working out the algebra, the field equation for ip a gives no new information. 
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6.1 SU(2) 

In SU(2) we have a pair of one-component weights. Taking S = e llfi ^ Tl , the ansatz 
for the Higgs fields, and the equations of motion are (h\\ = hi), 

ip 1 = h 1 T 1 , , ip ai = hST^S' 1 , , ip ai = hST^S' 1 , (130) 

^O^-a") =g{l-a)h 2 , (131) 
h" + (1/pK = p 2 h + {k + Xh^h 2 , (132) 

h" + {l/p)h! - (1 - a) 2 (l /p 2 ) h = p 2 h + Khh x + \{h 2 + hi) (133) 

This is essentially a generalized set of Nielsen-Olesen equations extended by the pres- 
ence of a profile function hi(p) that is not required to vanish at the origin, and tends 
to v c , when p — > oo. 



6.2 SU(3) 

For definiteness, let us consider a vortex with weig ht /?! = (y^/2,1). Using that, 
«i • /3i — 1, «2 • Pi — 0, and «3 • /?i = 1, the off-diagonal fields in eq. (122) become, 



V'ai = h ST^S , ?/'q, 2 



h T ( 



<>2 



(134) 
(135) 



The use of the same profile function for ip ai and tp a3 is necessary in order for the left 



and right hand sides in eq. (127) be oriented along the same direction in the algebra. 



In this case, the sum over a gives, 

2g (a ■ (3 q )(l — a)h 2 a a ■ T 



and replacing in eq. (127), we get, 

1 /V 
9 V P 



2g(\-a)h 2 (ai + a 3 )-f 
g{l-a)h 2 Px-f , 



g{l-a)h 2 . 



(136) 



(137) 



Now, when summing over positive roots, we have a| g a| p = hS qp , and the field 



equation (128), associated with Higgs fields along diagonal directions, turns out to be, 
d 2 h qp = p 2 h gp + nh 2 5 qp + Xh 2 h qp 



+2(na 2 \ q + Xa 2 \ q >h qql )(h 2 - h 2 ) a 2 



p ■ 



:i38i 
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The equation for ip a leads to, 

d 2 h a -(l-a) 2 {a-p 1 ) 2 {l/p 2 )K = 
= p 2 h a + 2Kh a h qq ia\ q a\ q ' + (A/3) h 3 a + (A/6) h a ft 



2 

7 

7^a 



hqqn (139) 

(with summation over q, q', and no sum over a). The a± component implies, 

d 2 h- (1 -a) 2 (l/p 2 )h = 
= p 2 h + 2Khhqq>di\qai\ q > + (A/3) h 3 + (A/6) h(h 2 + h 2 ) 

+ (k/3) hh + Xa 1 \ q iai\ q //hh qq/ h qq// . (140) 

while for the a 2 component, we obtain, 

d 2 h = p 2 h + 2Kh h qq >a 2 \ q a 2 \ q > + (A/3) hi + (A/3) ft ^ 2 

+ (/t/3)/i 2 + Xd 2 \q'd2\q"h h qq/ h qq /> . (141) 

The equation for a 3 is similar to that for c?i, with a 3 in the place of «x. 
Now, considering the following structure, 

h qp = -^hi j3x\ q /3x\p + 3h 2 a 2 \qd 2 \ p , (142) 

and using, 

= - + 3 2 |g«2|p , (143) 

the equations for the Higgs profiles accompanying frame components that have no 
defects become simplified to the form, 

hi + (l/p)h' = p 2 h + (re/3) (2h h 2 + h 2 ) + (A/3) h (h 2 + h 2 + h\) , 
h'[ + (l/p)h[ = p 2 h x + Kh 2 + \h 2 hx , 

h' 2 ' + [l/p)h' 2 = p 2 h 2 + Kh 2 + Xh 2 h 2 + (2/3) (k + \h 2 ) (hi - h 2 ) . (144) 



From eq. (143), the boundary conditions for h qp amount to hx —> v c , h 2 — > v c , when 
p — > oo. 

Finally, the equations for d\ and d?3, that correspond to the frame components 
having defects, become simplified to a single equation for the profile function h, 

h" + (l/p)ti - (1 - a) 2 (l/p 2 )h = p 2 h 

+(«/6)/i (2ho + 3hx + h 2 ) + (X/12)h (6 h 2 + 2h\ + 3/i 2 + h\) . (145) 

Therefore, we obtain again a modified set of Nielsen- Olesen equations, but coupled to 
three profile functions h ,hx,h 2 , which are not required to vanish at the origin. 
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7 Junctions 



In this section, we shall write the non Abelian phase factors corresponding to different 
types of junctions. They will provide the guiding centers, and the charges, or fluxes, 
of center vortices and monopoles. These fluxes will be defined with respect to the 
outward normal to the surface S 2 around the center of the junction. These factors will 
also define those gauge and Higgs field components that have defects. Therefore, they 
are fundamental to obtain the appropriate boundary conditions, as well as determining 
the type of quark sources they can bind. 

Before doing this, we would like to comment about the smoothness of the gauge 



field behavior in the asymptotic monopole region, given in eq. (91), and the form of 
the monopole density current. As the off-diagonal fields Cf depend on frame compo- 
nents containing defects, it is important to check whether or not these singularities are 



cancelled in the second term, C^n a . From eq. (49), and using 

R~ l diR\ B c = e B ■ R~ l diRe c = n B ■ dih c , 

we have, 

CfM A \ BC = -Remise , -gfABcCf = (n B ,d t n c } . (146) 



9 



Then, we can write, 



n q A difip = (n A , <9;n p ) n q An A = -g f DAp Cf f qAB n B , (147) 

and defining, 

Pp = $\*n q = S$.fS- x , (148) 

— * 

where the /3's are (2N times) weights of an irreducible representation, we obtain, 

PpAdiPp = g0-M)X h C?n b . (149) 
Now, for SU (N) we have the property, 

J20-M) 2 \ ab = c r 5 ab , (150) 

where c r only depends on the irreducible representation of SU(N) considered; for the 
fundamental representation, c r = 2. This can be verified by acting with both members 
on the nonzero root vectors, and using the Freudenthal - de Vries formula for J^^(w-a) 2 
[561 [57]. Thus, we get, 

C?n a = {l/gc r )Y,Pf>*diPf>, (151) 
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that is, the gauge field in the asymptotic monopole region (cf. eq. (91 )) can be written 



in terms of A\ and the diagonal components n q , which are regular on that region. Note 
also that the set of weights is invariant under the Weyl group (221 G3] • 

We can proceed in a similar manner to get information about the monopole density 



current h\ in eq. (96), that is given in terms of H\- in eq. (95). From the general 
property, 



(no, diU A A djU B ) = (nc, n A > A n B ')(n A >, din A )(n B >, djn B ) 
= ig 2 M A M c M B \ DE C D C* , 



(152) 



where we used the algebra and eq. (146), we obtain 



(Pp, A djPp) = ig 2 ■ Mf\ DE C D C E , 



fj, V 



(153) 



For a fundamental weight /3, we have, (/3 ■ M) 3 = (3 ■ M, as both members gi ve th e 
same value ■ a) 3 = f3 ■ a when applied on the nonzero root vectors (cf. eq. ( 116[ )), 
and vanish when applied on the zero root vectors. Therefore, we get, 



0|«fl5 = -gP\ q f qde CZCZ = — (Pp, diPp A djPp) . 

g 



(154) 



That is, the projections of the densities only depend on the coset, they are smooth 
everywhere but at the pointlike defects in the local diagonal sector, and have the same 
form as the well-known topological charge densities for S 2 — > S 2 mappings. 



7.1 iV-branch normal junction 

The normal junction is simply formed by assembling center vortices or confining strings. 
In £j8j we shall see the natural way to introduce quark probes in the dual model. After 
this, it will become clear that this type of junction confines a color singlet state of iV 
external fundamental quarks to form a normal baryon. The normal junction can be 
defined by the mapping, 

g — e *(vi 01+-+VN-1 Av-i)-? (155) 



where, for q = 1, . . . , N — 1, the weights (3 q are given by eq. ( 115 ), and the multivalued 
phases (p q satisfy, 

J q \i = £ijkdjd k ip q , J q = J ds^ §( 3 )( x - x q (s)) , (156) 

x q G -R 3 . For each path, the parameter s runs from — oo to +oo; in these limits the 
paths tend to infinity in R 3 , and they are chosen as follows. In the interval s G (— oo, 0] 
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all the paths coincide, describing the N-th. branch, running from infinity to the point 
O at s = (the center of the junction). A surface S 2 around O is intersected by this 
branch at a guiding center c^. Then, for parameters in the interval s G [0, +00), the 
N — 1 paths separate forming the other N — 1 branches that intersect S 2 at the guiding 
centers c q , q = 1, . . . N — 1. 

The mapping R = R(S) is ill-defined at the centers Ci, . . . , cjy. On the disks D Cq 
centered at c q , q = 1, . . . , N — 1, S(x) behaves as, 

S Cq (x) ~ S Cq eWtar* , S Cq = J] e***t*' f (157) 

where the factor S c is smooth, as it depends on phases ip p , p 7^ q, which are regular at 
c q . On the disk D CN centered at c/v, when its radius e gets smaller, all the phases have 
approximately the same value, tp q ~ ip, q = 1, . . . , N — 1, thus implying the following 
behavior, 

S CN (x) ~ e iv(0l+-+0N-i)-f = e -i V ?N-T _ ( 15 g) 



Then, comparing with eq. (79), we can identify, 

Xc,~¥>, , Pc q = P Q , q = l,...,N-l, (159) 

Xc N V 5 , /3 CJV = Pn ■ (160) 

The center vortex charges, labeled by {3 Cl , . . . , (3 CN , then correspond to the N weights of 
the fundamental representation, satisfying the condition /3 Cl + . . -+(3 CN = 0. This means 
that there is no monopole at the origin, and the junction will match N fundamental 
quarks having the N charges different. 

It is also worth mentioning that the same mapping S can be rewritten in terms of 
any set of iV — 1 independent fundamental weights, by simply changing variables to a 
set of phases whose associated currents coincide on a different branch. 



7.2 Two-branch vortex-monopole junction 

This junction is formed by a pair of center vortices attached to a monopole, one of the 
vortex charges belongs to the fundamental representation, while the other belongs to 
the antifundamental. When adding external quark sources in the dual YMH model, this 
junction will be able to confine a color nonsinglet quark/antiquark pair that becomes 
bound to a valence gluon to form a hybrid meson state. For general SU(N), it can be 
defined by the mapping, 

S = e*f* f W(x) , W(x) = e i6V " T "«p , (161) 
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where (p and 9 are the polar angles, and the roots and fundamental weights are given 



by eqs. (103), (115). At the north pole, the behavior is, 



S Cl (x) 



(162) 



At the south pole, W(x) approaches W a = e 



iir VNT a 



a Weyl reflection whose effect is, 



W~ X P'T W a = /3'-f. , (3' = (3-2a{a ■ (3)/a 2 

In our case, where (3 = (3 q and a = a qp , 

ft = ft q - 2Na qp (a qp ■ f} q ) 

= q - 2N(w q - w p ) = (3 P , 

so that the behavior at the south pole is, 



(163) 



Then, comparing (162) and (165) with eq. (79), we obtain, 

Pci Pq i Pc2 Pp i 



and the monopole charge is (cf. eq. (99)) 

2tt 



Q, 



9 



C2 / 



2tt 
9 



2N{w q - w p ) 



2tt 
9 



2Na 



qp ■ 



(164) 
(165) 
(166) 

(167) 



As the roots a qp generate the lattice A(Ad(SU(N))), the monopole in this juncion is 
expected to be an elementary object of the dual model. In fact, it can be naturally 
identified with a confined valence gluon with color a qp . 
In the local basis, 



Tlq — STqS 



n r 



ST^S- 1 



n. 



STsy S 



'1681 



the components having a vortex-like defect at 9 = (9 = it) are those n a , such that 
T a does not commute with {3 q ■ T • T). This is the case when one of the fundamental 
weights forming a (cf. eq. (103)) coincides with w q (w p ). 



In the case of SU(2), we can take for example, 

g _ e iyV2T le i8V2T ai 

which leads to, 

rii = cos 9 T\ + sin 9 (cos (p T &1 + sin ip T ai ) 
Uocx = — sin 9 T\ + cos 9 (cos <p T &1 + sin T ai ] 
n ai = -simp T Sl + cos y? T ai , 



(169) 



(170) 
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and the components n ai , n &1 rotate once when we approach the positive and negative 
z-axis. This is the frame depicted in Fig. 2. 
In the case of SU(3), taking, 



(171) 



the components n a ,n a with a = cti,a 3 (resp. a 2 ,<3 3 ) contain a vortex-like defect on 



the positive (resp. negative) z-axis. For the definition of ai,a 2 ,a 3 , see eq. (108). 



For general SU(N), the local diagonal directions can be analyzed by considering 



the projections Pp in eq. (1 



addition, because of eq. (101 
f3 q and j3 p , 



8). Among them, those with /3 ^ P q , /3 P are trivial. In 
I, (3 q + (3 p can be written in terms of weights other than 



so that the combination Pp 



+ Pp p is also trivial. Then, in the local Cartan 



sector, the nontrivial frame dependence is contained in the combination, 



S((3 q - f3 p ) ■ TS 



-i 



2L r 



a 



a 



(172) 



a 2 



(173) 



Note that the Lie algebra element L a , is a local version of the diagonal generator of 



the so (3) subalgebra in eq. (33). Using eq. (161), we obtain 



a " ' T i • a ( T « , • Ta 
cos V — - — hsinfi I cos ip + sin ip —= 



'Of 



(174) 



In other words, the frame defined by eq. (161) describes a center vortex with 
weight (3 P entering a surface S 2 around the origin, a center vortex with weight (3 q 
leaving, interpolated by a monopole configuration. The latter correlates S 2 , the space 
of directions around the origin, 



f = (cos 9, sin 9 cos ip, sin 9 sin (p) , 
with L a , a generator of a local so (3) algebra labeled by a = a qp , 



rif 



'of 



(175) 



Note also that at 9 = and 9 = tt, the ill-defined phase p is not manifested in the 
local Cartan sector, as required by the general discussion given in £|4j In that sector, 
the only defect is at the origin, where L a is ill-defined, so that the profile functions 
accompanying L a will be required to be zero at that point. 
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The frame dependent fields in eq. (43 ) are important for the general parametriza- 
tion of the gauge field in eq. (46). Using eq. (49), we get, 



CfM A 



2g 



di<p + cos 0(&-& 
(di<p sin 6 Ma -di9M a ) . 



■M 



(176) 



As expected, when 9 — > (r esp. 9 — > n ), the diagonal part is oriented along (3 q (resp. 
(3 P ). In addition, from eqs. (148), (151), we have, 



Cfn a = {l/2g)[P Pq A d^ q + P Pp A d r P Pp ] = (1/g) L a A d r L a . 
Therefore, we get, 



A t = (A?-Cf)n A 



- Cf n A 



Pfk 



diip 



[I + COS I 



+ Pi 



' 1 — COS I 



/3 P " 



1 

9 



L a A diL a 



(f3 q + (3 p )-T + L c 



1 



{(3 q + (3 p )-T-L c 



(177) 

(178) 
(179) 

(180) 



The second term in eq. (179) generalizes the one obtained in the asymptotic region 

This equation also puts in evidence the 



of an SU{2) 't Hooft-Polyakov monopole 
guiding centers of the junction, as well as the boundary conditions to be satisfied by 
Af. Besides the conditions Af — > in the asymptotic region, and A% — > in the 
asymptotic monopole region, these fields must compensate the singularities present at 
9 = and 9 = n, in the first term of eq. (179), and the singularity at the origin, in 
the second. 

For the projections of the monopole current Hf- in eq. (154), we note that those 
with (3 7^ (3 q , (3 P vanish, as in that case Pp is x-independent, while using eq. (180), 

1 



Mm 



9 



{Pf3 q ,diL a A djL a ) 



(PfrALahdjL*) . (181) 



Then, taking into account eq. (174) we also obtain, q + /3 p \ g ' Hfj = 0, so the nontrivial 
part of the monopole density current is contained in the projection along the root a qp , 



a\ q H q l3 /a 2 = --(L a , d { L a A d d L a ) 



a = a 



qp ' 



;i82) 



where we recognize the topological charge density for maps from S* 2 — > S 2 , namely, 
from the space of directions around a point into the space of local diagonal directions 
in a local so (3) subalgebra. 
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7.3 iV-branch vortex-monopole junctions 

These junctions will confine a general color state of N fundamental quarks, binding 
them to an appropriate monopole, forming a hybrid baryon. They can be defined by 
a mapping of the form, 

s = e <(<w N - 1 yf w ^ ^ ^ 183 ) 

where the first factor is the one describing an iV-branch normal junction. 
On the disks D Cq centered at c q , q — 1, . . . , N — 1, S(x) behaves as, 

S Cq (x)~S Cq e l ^P«- f W(c q ). (184) 
On the disk D CN centered at c^, when its radius e gets smaller, we have, 

S CN (x) ~ e ^(A+-+/3iv- 1 )-T W ( CN j = e ixc^r%( Cjv ) _ (185) 

Then in order to have smooth mappings along the local Cartan directions, n q = ST q S~ l 
on S 2 , we shall require H^(ci), . . . , W(cn) to be either the identity or a Weyl reflection. 
The monopole charge is proportional to the sum ^ Yl c Pc, which being the sum of N 
weights of the fundamental representation, a reflected subset plus an unreflected one, 
is necessarily in the weight lattice of Ad(SU(N)). 

8 Center vortices and chains in pure YM 

As we have seen, in the context of Yang-Mills-Higgs models, center vortices and junc- 
tions are interesting objects that could describe the states of the gluon field needed to 
confine quarks, and form the different normal and hybrid hadrons. To complete this 
picture, it is important having a description of heavy quark probes in the language of 
the proposed YMH model. This leads to another question, namely, why there could 
be some relationship between the YMH model and the pure QCD Yang-Mills theory 
it is supposed to describe. Of course, this type of question corresponds to the open 
nonperturbative problem of confinement, so that an answer from first principles is out 
of reach. What we could do is, based on some confinement scenarios in the lattice, 
suggest some natural phenomenological choices, once the dual model is assumed to 
have the form given in eqs. ([I]), ( [18] ). 

In the introduction, the description of the thin configurations in SU(2) YM theory 
motivated the analisis of topological smooth configurations in the YMH model. Now, 
we can use the general properties of the latter configurations, discussed for a gauge 
grop G, to write those of the former, for a gauge group G c . Moreover, we would like to 
suggest, by means of a heuristic argument, a possible relationship among both sides. 
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Let us consider a YM theory with gauge group G c and coupling constant g e , in AD 
Euclidean spacetime, 



S e YM = J d 4 x £ -(F; u ,F^), (186) 

for a gauge field A e ^ representing gluons. As is well-known, its analysis in the lat- 
tice supports scenarios dealing with center vortices, monopoles, and chains formed by 
them. For example, some observed properties of the quark potential, such as asymp- 
totic linearity and iV-alitjQ are well described by the ensemble of center projected 
configurations (see refs. [2j [3] ) . They are obtained by initially choosing a gauge such 
that the link variables U^x) are as close as possible to some Z(N) lattice configuration, 
then decomposing, 

U^x) = V^x) Z^x) , (187) 

where Z^[x)l is the center element of SU(N) closest to U^x), and finally discarding 
the "perturbative part" V^x), 

U^x) -»■ Z^x) . 

The obtained configurations essentially correspond to center vortex defects, located 
at plaquettes where J^Z^x) ^ 1. At low temperatures, center vortices percolate, 
and the average of the Wilson loop in the ensemble displays an area law, that is, 
confinement. While this scenario is very good at describing iV-ality, the stringlike 
behavior of the confining potential, that is, a nontrivial Luscher term associated with 
transverse fluctuations of the string, has not been observed in the center-projected 
data. 

In the lattice, monopoles and antimonopoles are also detected. Most of them have 
been found to be correlated with center vortex worldsheets to form chains, configura- 
tions that are considered very promising candidates to explain the different properties 
of the confining quark potential. Scenarios only involving an ensemble of monopoles 
have also proven successful in a number of respects. In this case, the confining phase 
corresponds to a monopole condensate. For a detailed discussion of these topics, see 
(2j 13] and references therein. 

Similar considerations to those given in £|3] and ^4] for the YMH model, here will 
lead to center vortices concentrated on worldsheets, that can concatenate looplike 
monopoles to form vortex-monopole thin junctions, or chains, in the language used in 
the lattice literature. Besides the dimensionality, there are important differences with 
respect to the YMH configurations. At the classical level, the YM theory has no scales, 
at the quantum level, smeared topological defects will be postulated to be relevant field 



configurations to compute averages. In this regard, a configuration of the form (46), 



3 the fact that in SU(N) YM theories the asymptotic string tension depends on the behavior of 
the quark representation under Z(N). 
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with — > "Pf, ua — > ua, representing a perturbative and nonperturbative sector, 
is the analogous of the lattice decomposition (187), generalized to include all types of 
defect, 

Al = {Vi-Z*)u A , Z* = -(l/g e )r ABC (u B ,d lUc } , u A = ST c A S-\ (188) 

where S G G e and the T^'s are generators of the Lie algebra $j e - The field strength 
tensor is (cf. eq. (50)), 

F; v = G$ u A , Git = (F*(V) - F*{Z)) . (189) 

Here, it is understood that quantities depending on the fields and Z^, involve the 
chromoelectric structure constants f ABC and coupling g e . 

For example, a closed center vortex worldsheet £ is described by, 

S = e*** , ^veA(Gl), (190) 

where x is a multivalued function, changing by 2tt when we go around a path linking 
S. Using eq. (73), we have, 

-F% f (Z) = -fj\"[d lt ,d t ] X , (191) 

Qe 

and the vortex guiding centers are represented by a contribution to the dual tensor 

= — rfr l^ v 5^\x-y{a u a 2 )), (192) 
d = da 1 da 2 tt-^t; ^^T~ > ( 193 ) 



where d 2 a fJiU integrates over S, parametrized by 27(<7i,er 2 ) (see |H]-[l9]). 

For vortex-monopole junctions, the mapping S has the general properties discussed 
in §|4j so that, in the dual tensor, there is a sum over open center vortices attached to 
looplike monopole worldlines. To identify the monopole guiding centers, note that for 
each open vortex, we have, 



2tt f 

-d u ^ u (Z)\ vort = —rj\« d 2 a^dj^(x-y(a u a 2 )) 

9e J 



yV\ 9 (j> dy^ <5 (4) (x-y)- j> dy^ 5 {A) (x - y) 
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where C + (C ) is the loop where the monopole (antimonopole) is localized. On the 
other hand, at say the monopole, the weights of the attached vortices must satisfy 



a condition analogous to eq. (83). That is, for each monopole (antimonopole), the 
guiding centers are represented by a nonzero contribution to the divergence of the dual 
tensor, 

- d u ^ u (Z)\ mon = ±— (ff Cl + ff C2 + . . .)|« / dy, S<Q(x - y) , (194) 

^ (ff cl + ff C2 + . . .) e A(Ad(Gl)) . (195) 

In YM theory, when the vortex guiding centers leaving a monopole are superim- 
posed, they form an unobservable open Dirac worldsheet. In this case, the only objects 
expected in the ensemble are the monopole worldlines. 

8.1 Mesonic Wilson loop 

For a given irreducible representation of G e , the Wilson loop for external charges is, 

W c [A e ] = (l/d r )trPexpU# e jf^AM 

= (l/d r )tr[Sf P exp ^ig e j> dx^ V^j S, 

= i(C)W c [V], (196) 

where Si and Sf are the initial and final values of S when we go around C, a loop 
that represents a process where heavy quark antiquark probes are created, propagated 
and then annhilated. Here, we used that the product S^ 1 Sf must be an element 3(C)/ 
in the center of G e , depending on C, the distribution of defects, and the irreducible 
representation used for the external quarks. Regarding this factor, let us initially 



consider it in the presence of a closed vortex defined by eq. (190). When the center 
vortex is linked, 3(C) gives e t2nr, ' Tc G Z(G C ). This center element is a number times 
the identity matrix, so it can be computed as the eigenvalue obtained when operating 
on any weight vector of the quark representation, labeled by a weight w e . That is, 

where L(C, E) is the linking number between C and S. In addition, the linking number 
can be equated to the intersection number I(S(C),T,) between S(C), a surface whose 
border is C, and the surface S, 

I{S{C),E) = l -j d 2 ~a, u j d 2 a, v 5^\w(s,T)-y(o u o 2 )) , (198) 
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d *V = -^drds — - — — j , (199) 

where w(s,t) is a parametrization of S(C) JIT} US]. Then, introducing the source, 

d 2 a^<5 (4) (x- w(s,r)) , (200) 

5(C) 



and using eq. (192), we obtain, 

= e -|/d 4 ^9oS, 11 ,Wo-/ M1 ,(Z) _ (201) 

This representation can also be used for vortex-monopole junctions. In this case, 
there is also a (richer) concept of linking, as they are formed by monopoles and an- 
timonopoles concatenated by center vortices to form closed objects. In addition, the 
intersection number can be applied to open surfaces, i.e., S(C) and the open vortices 



in the complex. Then, the center element 3(C) can also be computed using eq. (201), 
that because of the flux matching between monopoles and center vortices, will give the 
same result for any possible S(C), as long as the monopoles do not touch S(C). 

8.2 SU(N) baryonic Wilson loop 

Here, we shall consider G c = SU (N), and the variable for a state formed by N funda- 
mental quarks. The baryonic potential can be measured from the average of the gauge 
invariant baryonic Wilson loop [1]-[E], [55] . 

W Nq [A e ] = ^ ,e ai ... aN e bl ... bN U aibl (A c ^i)--- U aNbN (A e ,T N ) , (202) 
in a similar manner to the qq potential, that can be measured from the large distance 



behavior of the average (Wc)- The U- factors in eq. (202) are the matrix elements of 
the holonomy, 

U(A e ,T) = Pexp (ig e [ dx.A e J , (203) 



\ Jr J 

and the paths r 1; . . . , represent N worldlines starting at a common point x, where 
the iV quarks are created, to a common final point y, where they are annhilated. 

Using the composition property of the holonomy, we can compose the different 
quark worldlines with a common path Tq 1 , that goes from y to x, thus forming iV 
closed paths C\ — T% o T 1 . . . Cn = IV 1 , and write, 

U(A c ,T)\ ab = U-'iA^r^^UiA^r oT^lcb . (204) 
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Now, for the different worldlines, the first factor corresponds to elements of the same 
matrix, so replacing in eq. (202), and using det U~ x = 1, we get, 



^W^ 6 ] — -JTj £a 1 ...a N £b 1 ...b N U aibl (A e ,Ci) . . . U aN b N (A e ,C N ) 



(205) 



In addition, for each closed path, we have, 

S,(C)= i {C)S i , 
U(A C , C) = S f (C) U(V, C) S- 1 = 3 (C) S t U(V, C) . 



Then, replacing in eq. (205), we obtain, 



W Ng [A e ]= i (C 1 )... i (C N )W Nq [P] 



(206) 
(207) 

(208) 



where 3(C) is given by eq. (201 ), using as w e a weight of the fundamental representation 
of SU(N). This formula generalizes the one obtained in ref. [58] for an ensemble of 
closed center vortices. 



8.3 Natural identifications 



The use of the decomposition (188), (189) and the observables (196), (208), pose a 



difficulty, when setting the fluctuations to zero, the thin defects have divergent action. 
On the other hand, the success of lattice scenarios involving defects could indicate 
that the proper ensemble is in fact formed by thick magnetic objects, with dimensional 
parameters depending on Aqcd- However, determining if this is the situation or not 
constitutes a difficult nonperturbative problem. At the perturbative level, the one-loop 
calculation of the effective action for a homogeneous magnetic diagonal background is 
lower than zero [59], so magnetic configurations seem to be favored with respect to the 
classical vacuum. However, there is an imaginary part in the calculation [00], implying 
a decay into off-diagonal gluons. This instability also occurs for vortex backgrounds 
[6T1 [50], and it seems to be a general feature of magnetic bakgrounds. In ref. [52], we 
pointed out a possible way out for the instability problem, reexamining how to define 
the coupling between a center vortex and quantum fluctuations. The idea is that 
fluctuations should not simply "see" the vortex as a diagonal background, but they 



should be coupled by the parametrization (188), where the presence of frame defects 



in the off-diagonal sector would change the counting of bound states of the fluctuation 
operator. 

The difficulties associated with a formulation from first principles of scenarios in- 
volving magnetic defects in the continuum led to different models to compute averages, 
after postulating some possible properties of the ensembles. For example, an effective 



40 



model of random surfaces, describing an ensemble of closed thick center vortex world- 
sheets, has been proposed in ref. [17]. These procedure implements the lattice idea of 
center projection. Arguments in favor of a Y-shaped potential between quarks were 
given in ref. [8], assuming closed center vortices, and the statistical independence of the 



total linking numbers through C\ . . . Cn in eq. (205). In ref. |17j . a dual Abelian YMH 
model for 577(3) YM has been proposed. It contains a Higgs sector that provides 
an effective description of an ensemble of magnetic monopoles. This type of model 
has been derived by assuming Abelian dominance, following a dualization procedure 
similar to that given in refs. [HI [63] for compact QED(3) and QED(A) (see also [62J). 

A general ensemble should include the possibility of vortex-monopole junctions. 
However, it is no clear how to obtain an effective field theory when the defects are not 
one-dimensional. This is the case of junctions in 3D Euclidean models [53J. In 4D, 
center vortices are two-dimensional worldsheets, for this reason, in order to suggest the 
natural choice of external sources and gauge group in the non Abelian YMH model, we 
shall only discuss the effect of monopoles, attached to unobservable Dirac worldsheets. 

Let us sketch the main steps that can be followed after assuming Abelian domi- 
nance, and including mesonic and baryonic variables. Disregarding off-diagonal fluc- 
tuations in Vfj,, the holonomy does not require path ordering. Then, the eigenvalues 
needed to compute traces and determinants involving U(V, C), can be simply obtained 
by acting on the weight vectors of the quark representation. In particular, the mesonic 
Wilson loop results, 

W C [A C ] = (l/d r )J2 e ^ fd * X9eS ^™ e<f '*" i ' P) ~ ff>v{Z)) , ( 209 ) 

where the sum is done over the weights of Q e , corresponding to the irreducible repre- 
sentation considered. For a baryonic SU(N) state formed by N fundamental quarks, 
we get, 

WjvJA 6 ] = — y~]e?f dix 9^l v w P1 +...+s^ v iS PN )-{fi^ v {v)-P llv {z)) ^ £210) 
p 

where the sum is done over the N\ permutations of the weights of the fundamental 
representation. 

The different terms in these expressions are expected to have the same average, 
after integrating over the ensemble with an appropriate phenomenological measure, 
and over the gauge fields, using the weight, 

£ - f d*x %3% U -G% _ e -J d±xl(F^(V)-F^(Z)f _ ^ll) 

Linearizing the kinetic term, both averages have the form, 
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J, 



/1U ) 



pV 



ge/3 e s 

gJi s \u + ■■■ + gJn sJL , 



for 
for 



W 
W 



W c 
W Nq , 



with ( = 1/AN e. In eq. (212), the integration variable A M „, which is coupled to dual 



field strength tensors, has been rescaled so as to work with external currents having 
fluxes with normalization g e (3 = 2Nw. To simplify the discussion, we shall initially 
assume that monopoles are singular pointlike objects. Note that the variable is not 
required to satisfy any boundary condition on the unobservable Dirac worldsheets, so 
we can directly integrate to obtain the constraint, 



^pvpa&v^-po 

Its solution leads to the dual representation, 



(213) 



(W) 



(214) 



where Gj„ 



d u A^, and the index P refers to fields in the Abelian projected 



case. The ensemble average is not well defined unless some phenomenological dimen- 
sional parameters in the integration measure be introduced. Different choices for the 
string tension and interactions among the monopoles will lead to a variety of effec- 
tive field models, exhibiting different phases |M]-[66]. Here, our main concern is only 
to suggest some natural identifications, which can be even seen after turning off the 
interactions. 

Using eq. (194), including the sum over the different monopole worldlines, we see 



that the dual vector field is minimally coupled to particles whose charges are given 
by (2n/g c ) times an element of the lattice A(Ad(GD). Considering the elementary 



monopoles in eqs. (194), (195), which correspond to the roots of Q v e , we have, 



(W) = / [DA P ] e 



7 7 



(215) 



1 



N 



e ah n e 

N ' k=l 



a- A 



(216) 



Here, we have set the interaction among the monopoles to zero, so that Z a , the 
partition function for monopole loops with charge a e v , is obtained by summing over 
the number of loops, and for each sector, path integrating over all possible loop shapes 
with the measure [-Dx^]. To avoid overcounting, only the positive roots need to be 
considered. The configurations are weighted by the interaction in eq. (214), between 



the monopole loops and the projected dual field, as well as the cost to have a loop with 
length L, depending on a phenomenological parameter m. 



42 



Now, as shown in refs. |67J El], the right hand side in eq. (216) is a functional 



determinant for a vacuum to vacuum amplitud in "particle" representation. This 
means, that Z a can be cast in the form, 

Z a = J[D<t> a )[Dt a )e-f^- D * +m2 ^ , D« = d»-i— a v e -A*. (217) 



where <p a is a complex field associated with a root « e v . In this respect, including density- 
density monopole interactions, some effective models have been obtained where a SSB 
phase is observed for m 2 < |62j. 

Summarizing, when considering Abelian projected configurations, and after turning 
off the interactions, the Wilson loop is, 

(W) = j [DAr][D<i> a ][D$ a ] e -^ rf4a;[ 4^-^) 2 +E a [l^^l 2 W^^] _ ( 21g ) 

In the case of the mesonic Wilson loop, that has been written for a general g e , a e v sums 
over the positive roots of q v c . For the SU(N) baryonic Wilson loop, the sum is over 
the N(N - l)/2 positive roots of su(N). 

Then, for the class of models we presented in £j2j it is natural to introduce the quark 
sources in the form, 

S = J d A x (^(Fpv ~ Jf+u, F^ v - J^y) + ^(D^ipA, D^4> A ) + Vm sg8 (ipA)j , 

and require that this action, when restricted to Abelian configurations, coincide with 
that derived for the Abelian projected effective model. Considering, 

A^ = Al-f , J^ = J^\ q T q (219) 

— -j= — = 4> a E a , — -j= — = (j) a E— a , (220) 
the covariant derivatives in 5" depend on the roots of g and the coupling g, thanks to 



the fact the Higgs fields are in the adjoint, while in eq. (218), they depend on the roots 
of e v and the coupling 2n/g e . In addition, the effective Abelian model in eq. (218) 
contains as many complex fields as positive roots of the Lie algebra. Then, from this 
perspective, it is natural identifying, 

2tt 

= 0^ , 9= — , (221) 

9e 

so when the interactions are turned off, the sector of off-diagonal flavors in the Abelian 



projected action S reduces to the exponent in eq. (218), with pi — m 
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8.4 Induced topological configurations 

Note that the parameter ( does not change the values of the center vortex and monopole 
charges of the YMH model discussed in ^3] and j |4} They are still given by eqs. (67), 



(71), and (f77|, so after the identification (221), that implies G v = G c , the vortex 



charges become electric charges, 

Qv = gJ , ^^A(Ge). (222) 
These are precisely the fluxes that can match those of the external quark sources in 



eq. (8.3). 



In particular, in a process where a qq pair, or a q ... q baryon, is created at 73 in the 
far past, the charges are static at intermediate times, and then they are annhilated at 
Tf, in the far future, the parametrization of S(C) can be done as follows, 

W(T, S) = (t,x(s)) , Tj < T < T f , x e R 3 , (223) 

where x(s) joins the quark and the antiquark. This leads to static sources, 

soi(x) = , Sij (x) = -e ijk Jds^- 5 (3 \x - x(s)) , x G R 3 . (224) 

Changing the shape of x(s) amounts to changing that of S(C), an operation that 
should have no observable effects. In other words, x(s) can be thought of as a Dirac 
string where only the endpoints (the quark locations) are observable. If desired, for 
the mesonic Wilson loop we can use a pair of Dirac strings, one running from — zq to 
—00, the other running from +00 to +zq. That is, — K%j represents a flux (3 e leaving 
(entering) the quark (antiquark) position at — z (+z ). For the baryonic Wilson loop, 
as Pi + . . . Pn — 0, we can choose the N Dirac strings running from the quark positions 
to 00, so that —Kij represents fluxes Pi, ... , P N , leaving each one of the quarks. 

The energy functional for (dual) magnetic configurations in the presence of static 
quarks is 

E = Jd 4 x (^(fy - J ih Fij - Jij) + ^(D^a, Diip A ) + VkspWa)) , (225) 



and its minimization leads to eq. (52), with the substitution 



Fij -> ((Fij - Jij) , (226) 



while eq. (|53|) is unaltered. Now, we are interested in solutions that have smooth 

Jij. Writing and ipA in terms of Af, 



energy density, in spite of the singular term Ju. Writing A„ and ipA in terms of Af 
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riA, and Kab, as in eqs. (46), (54), and taking into account eq. (50), we then require 
a regular, 

(227) 



Jy = FMA) n A - F£(C) n A - J tj . 



In the mesonic case, this requirement can be fulfilled by inducing a finite length 
vortex with guiding centers x c (s), running from the quark to the antiquark, and weight 
(3 equal to the quark's weight (3 e . That is, by inducing a configuration parametrized 
by, 



S 



R 



(228) 



where \ is a multivalued function that changes by 2tt, when the closed curve formed 
by the composition of x c (s) and the string x~ 1 (s) is linked. Using the relation between 
couplings in eq. (221), the last two terms in eq. ( 1227) become, 



[l/g)[d i ,d j ]xPe-M + 



dlT^ — * — * 

ds=-8( 3 \x-x(s))p e -M = 

ds^5 {3 \x-x c {s))(J e -M . 
ds 



(229) 



This sum vanishes at the Dirac string x(s), so A^ must be smooth there and, as a 
consequence, no special condition for the profile functions h A B is required at x(s). On 
the other hand, at the vortex guiding centers, vanishes. Then, when x — > x c (s), 
the diagonal components of along /3 e must tend to —- g diip (3 e ■ T, for the sum of the 



first two terms in eq. (227) be smooth, In addition, when x — > x c (s), the Higgs profiles 



coupled to these components must tend to zero. These boundary conditions must be 
interpolated with those at oo, so the energy density will be distributed around x c (s). 

For the baryonic Wilson loop, we have seen that the external sources are associated 
with the N different weights pi, ... , (3^ of the fundamental representation of SU(N). 
In this case, a regular energy density can be obtained by inducing a configuration with 
Y-shaped finite guiding centers that, combined with the external sources, are described 



by the map S given in eq. (155). Again, because of the external sources, the usual 



regularity conditions only need to be imposed at the guiding centers, thus localizing 
the energy density on a finite Y-j unction. 

Then, when introducing external quarks, the equations of motion are those given 



in (52), (53) but with modified boundary conditions. A similar situation applies to 



the SU(N) hybrid meson. It can be induced by only using the regularity conditions, 
discussed in £|7.2l inside a sphere S 2 around the origin, and using a regular A^ out- 



side. The apparent singularities in the energy density outside S 2 will cancel against 
appropriate external sources, corresponding to a quark and antiquark with charges (3 q , 
-f3 p placed at — z and +z , respectively. 
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The search for the minimum energy must also consider a minimization with respect 
to the guiding centers. For example, for fundamental quarks in SU(3), when consid- 
ering a qq pair or three quarks on the corners of an equilateral triangle, this step is 
expected to be unnecessary, if the guiding centers form a straight segment joining the 
quarks, and the Y-junction branches are oriented at angles of 120°, respectively. 

Summarizing, the topological features of the different configurations, the boundary 
conditions needed to have a regular energy density, as well as the dimensional scales 
associated with the diagonal and off-diagonal Higgs sectors, are expected to lead to 
finite energy solutions to the equations of motion. However, it is important to empha- 
size that care must be taken so as to avoid spurious solutions. For example, in SU(N), 
following the steps we did for a qq pair, using a weight (3 e in the adjoint, it seems that a 
string confining adjoint quarks could be formed. The point is that we should not insist 
on screening the external sources only using the Abelian phase in eq. (228). When 
imposing the boundary conditions, all types of R mappings should be considered. A 



similar situation was discussed at the end of £ 3. 1 In fact, when the sources are in the 



adjoint, a non Abelian phase that screens the external sources, and is smooth on the 
region between the charges, will be induced, thus implying a vanishing string tension 
for adjoint quarks, in accordance with iV-ality. 



9 Conclusions 

In this article, we analyzed a class of Yang-Mills models containing adjoint Higgs 
fields, with SU(N) symmetry spontaneously broken down to Z(N), and showed they 
contain a rich variety of topological objects. These include center vortices, iV-branch 
junctions exclusively formed by them, and junctions where different center vortices are 
interpolated by monopole-like configurations. In the context of dual superconductor 
models, the first one represents the string that confines a color neutral quark/antiquark 
pair, to form a meson. The second possibility represents the generalized Y-shaped 
potential that confines a color neutral iV-quark state, to form a baryon. The third, 
that contains correlated monopoles, is a new kind of configuration that naturally arises 
in these models. The simplest can be thought of as a gluon field state containing a 
colored valence gluon represented by the monopole. They could confine colored quark 
states to form an overall color neutral hybrid object. For example, a junction formed 
by a pair of center vortices with different charges, interpolated by a monopole, can 
confine a nonsinglet quark/antiquark pair, to form a hybrid color neutral qgq' meson. 
Similarly, hybrid baryon states can be formed by iV-branch vortex-monopole junctions. 

These states are indeed allowed by QCD, existing some candidates in the observed 
hadron spectrum, lattice evidence in favor of them, as well as experimental collabora- 
tions searching for this type of gluonic excitation. 
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Although the main interest is on su(3), we have found it convenient to settle the 
discussion in the framework of a connected compact simple Lie algebra q, and use its 
general structure. In our model, vacuum configurations in the SSB phase are labeled 
by tia, a local basis for $j, that admits a Cartan decomposition into diagonal (n q ) 
and off-diagonal (n a ) local sectors. Equivalently, they can be labeled by a mapping 
R E Ad(G). 

As is well known, when the center of the universal covering group G is nontrivial, 
stringlike center vortices can be formed, as the first homotopy group of the space 
of vacua is given by Il 1 (Ad(G)) = Z(G). On the other hand, isolated monopole- 
like objects cannot, as U 2 (Ad(G)) = 0. Here, we have emphazised the fact that 
mappings from S 2 , the sphere around a point, into the space of vacua, only measures 
the possibility of having isolated objects. That is, the trivial IT 2 of a compact group 
does not forbid monopole-like objects attached to center vortices to form junctions. In 
such a configuration, on the sphere around the monopole, there would be regions of 
false vacuum, namely, the center vortex cores. As we have shown, these junctions can 
be formed indeed, with the monopole contained in the local diagonal sector of the Lie 
algebra. They are described by the exact homotopy sequence, 

= U 2 (Ad(G)) ->■ n 2 (D) ->• III (Ad (#)) -»• III (Ad (G)) = Z(G) , 

where D = Ad(G)/Ad(H) , and H is the subgroup of G leaving the diagonal generators 
T q invariant. In this sequence, the manifold D can be identified with the local diagonal 
sector n q = STgS -1 , S G G. The manifold H is associated with the center vortices, 
as when we get closer to their guiding centers they are essentially contained in a local 
Abelian subgroup. The S* 1 in Hi(H) refers to the composition of loops around the 
different center vortex cores crossing the surface S 2 , where n 2 (-D) is being analyzed. 
Then, when this composition is required to be trivial in Ad(G), an identification be- 
tween 112(D) and LTi (if) is established. This means that configurations exist where a 
monopole along the diagonal sector is attached to center vortices, as long as the charges 
of the vortices, which are essentially weights of the universal covering of the dual group 
G v , add up to a weight of Ad(G v ). In turn, the weights of Ad(G v ) become the possible 
monopole charges, and the absence of isolated monopoles can be interpreted as the 
confinement of valence gluons. 

In the case of SU(N), the simplest center vortices are known to be labeled by 
weights of the fundamental representation. On the other hand, the simplest monopoles 
are labeled by the roots, which can be written as the difference of a pair of fundamental 
weights, those associated with the pair of interpolated center vortices. 

As a possible model supporting the different topological objects, we proposed one 
containing as many adjoint Higgs fields ipA as the dimension of the Lie algebra, and 
constructed a potential displaying flavor symmetry with global group Ad(G). Its terms 
were obtained by considering Higgs field products in the Lie algebra and forming 
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invariants up to quartic order in the fields. When the cubic term is nonzero, as the 
square mass term is lowered, a first order transition to a SSB phase is obtained. When 
the parameter of the cubic term is nullified, there is a second order line. 

In the SSB phase, the minimization of the potential turns the set of fields ipA into a 
local basis ua, so that the previous analysis concerning the possible topological objects 
follows. The mapping R e Ad(G), associated with n^, is a non Abelian phase whose 
defects determine the guiding centers of the smooth topological objects, as well as the 
boundary conditions to be imposed on the gauge and Higgs fields. 

For SU(2) and 577(3) center vortices, we explicitly showed how to close the field 
equations for a reduced set of profile functions. The equation for a(p), the gauge 
field profile along the center vortex weight, is the same encountered for the Nielsen- 
Olesen vortex. Here, the Higgs profile h(p) is the one associated with those off-diagonal 
components that are rotated by diagonal transformations, with axis along the vortex 
weight, so it must vanish at the origin. The field equation for h(p) is a modified 
Nielsen-Olesen equation, containing couplings to profile functions for the Higgs fields 
that are not rotated. 

For general SU(N), we wrote the non Abelian phase describing the guiding centers 
of a junction formed by a pair of fundamental center vortices, with weights w q and w p , 
attached to a monopole. In this case, we showed how the monopole is contained in 
L a = a\ g n q /a 2 , a = w q — w p , the generator of a local 50 (3) algebra. For the gauge 
field, we put in evidence the presence of a center vortex in the positive (resp. negative) 
z-axis with weight w q (resp. w p ), and the monopole term L a A diL a . We also showed 
that the monopole charge is the topological charge for the mapping S 2 — > {L a }. 

The diagonal Higgs fields are important to set a scale for the monopoles, giving a 
smooth transition from the asymptotic monopole region, where the diagonal sector does 
not contribute, to the center of the monopole, where the field profiles accompanying 
the singular L a must tend to zero. 

The general picture is completed with a heuristic discussion about what could be 
the relation between the YMH model with gauge group G and pure YM with gauge 
group G c . The natural choice of external sources in the YMH model, to describe 
quarks in different irreducible representations of G c , is also addressed. The idea is 
that configurations involving monopoles, center vortices and vortex-monopole junctions 
(chains, in lattice language) also appear in pure YM theories defined on the lattice. 
The Abelian projection, leading to ensembles of monopoles, and the center projection, 
where only the center vortex degrees are maintained in the averages, proved to be very 
successful to describe some of the properties of the confining force. Chains are also 
candidates to explain the different aspects of this force. 

In continuum Euclidean YM theory, the different projections can be simulated by 
considering ensembles of magnetic defects, characterized by a phenomenological mea- 
sure that contains dimensional parameters. In AD, the monopole component is formed 
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by looplike defects, placed at the border of observable center vortex worldsheets, and 
unobservable Dirac worldsheets. As is well known, in Abelian projection scenarios, an 
effective description of this component is expected to be a field theory. This comes 
about as the monopole ensemble contains a sum over different numbers of closed world- 
lines, that is, a second quantized field theory. As the possible elementary charges of 
YM monopoles are roots spanning the lattice Ad(Gl), the Abelian projected effective 
model contains as many charged fields as the number of positive roots. These modes 
were identified with the off-diagonal flavors of the Abelian projected YMH model, 
when defined with adjoint Higgs fields transforming in Ad(G y c ), that is, it is natural 
identifying G = G y c . 

In this regard, it would be interesting to look for the ensemble whose large distance 
behavior is described by the complete YMH model we proposed, although obtaining 
an effective model outside the Abelian projection assumption is a difficult task, as well 
as dealing with ensembles of vortex-monopole junctions in AD. 

Because of the identification G = G v e , the YMH center vortex charges become 
weights of the universal covering of G v = G c . In addition, in dual language, we have 
shown that the mesonic and baryonic YM Wilson loops become naturally represented 
by external pointlike monopole sources, with weights of the irreducible representation 
of G c where the quark probes are defined. That is, the quark sources can be matched 
and confined by center vortices, junctions formed by them, and junctions containing 
interpolating monopoles. In the YMH model, these configurations are smooth, in con- 
trast to what happens in the YM theory. They are solutions to the classical equations 
of motion, representing normal and hybrid mesonic and baryonic states. Further re- 
search on this direction, performing a numerical analysis of the equations, will permit 
to obtain the model's potential between a quark/antiquark pair, for various irreducible 
representations, and the field distribution of confining strings and valence gluons inside 
hadrons. 
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